Journal of Molecular Biology

The SKP1/Cul1/F-Box (SCF) complex is an E3 ubiquitin ligase responsible for targeting a range of proteins for degradation by the 26S proteosome. Within this complex, a variety of F-box proteins (FBPs) link to the SCF complex via the SKP1 adaptor protein allowing for differential substrate recognition. In the intracellular parasite Toxoplasma gondii, SKP1 is subject to oxygen dependent regulation. Under normoxic conditions, the prolyl hydroxylase PHYa hydroxylates SKP1 priming it for modification by five SKP1-specific glycosyltransferase activities. Glycosylation plays an important role in controlling SKP1 by weakening the tight SKP1 homodimer and affecting the profile of bound FBPs in cells. However, the presence of the terminal SKP1 glycosyltransferase, GAT1, in the SKP1 interactome regardless of its glycosylation status is atypical for an enzyme. Furthermore, gat1-knockout cells exhibit a unique repertoire of FBPs bound to SKP1 relative to normal and other glycosylation-defective mutants. Utilizing sedimentation velocity analytical ultracentrifugation, we demonstrate that the native GAT1 homodimer complexes with SKP1 monomers with affinity and stoichiometry dictated by its glycostate. Computational modeling validated by mutational probing shows that GAT1 competes with the same core hydrophobic interface utilized by FBPs and the SKP1 homodimer. This interface is complemented by varying, transient fuzzy-like interactions contributed by the intrinsically disordered C-terminal region (CTR) of SKP1 that are in turn constrained by the glycan. Furthermore, substoichiometric levels of GAT1 mediate monomerization of SKP1 in a CTR-dependent manner, indicating that GAT1 has the kinetic potential to promote SKP1 monomer availability, with consequences on its FBP-binding preference in cells.

The Mitochondrial contact site and Cristae Organizing System (MICOS) complex is an inner mitochondrial membrane (IMM) assembly present at the cristae junction. It is responsible for regulating cristae formation and remodeling. However, its structure is not known. We applied Bayesian integrative structure determination to characterize the structure of the Mic60, Mic19, Mic10, and Mic13-containing MICOS complex combining AlphaFold predictions with data from crosslinking mass spectrometry, biochemical assays, electron tomography, homology modeling, and sequence alignments. The integrative structure revealed novel mutual interfaces among Mic10N,C, Mic60LBS1,LBS2,mitofilin, and Mic13central,C, which were experimentally validated. Several likely-pathogenic missense mutations also localize to these novel interfaces, highlighting their importance. Our results indicate that Mic13 likely facilitates MICOS assembly by binding Mic10 in the IMM-proximal region and Mic60 in the intermembrane space. Taken together, our integrative approach sheds light on the structure and assembly of the MICOS complex.

Short, linear sequence motifs within the human prion protein (PrP) may encode local aggregation tendencies that are not apparent from full-length sequence analysis. To map intrinsic amyloidogenic potential across PrP, we generated a complete one-residue-step library of overlapping 15-mer peptides from the 253-residue human PrP sequence and evaluated each peptide using WALTZ in both high-specificity and best-overall-performance modes, with full-length SNPeffect4/WALTZ output used for comparison. Peptide-level WALTZ analysis identified several candidate amyloidogenic regions, including an N-terminal signal-peptide segment spanning approximately residues 8-21/22 that was not detected in the publicly available full-length WALTZ/SNPeffect4 output. Predictions made using AlphaFold 3.0 also identified a short C-terminal area whose representative peptides formed either -helix or pair of {beta}-strands, indicating a region of potential higher susceptibility to conformational dynamics. This computational study supports the use of peptide tiling as a complementary screening strategy for identifying candidate short aggregation-prone motifs in PrP and other misfolding-associated proteins.

The androgen receptor is a hormone-dependent transcription factor that regulates a wide range of physiological processes and plays a pivotal role in the development of prostate cancer. Its 555-residue, intrinsically disordered, N-terminal domain is involved in the modulation of transcriptional activity by recruiting co-regulators and mediating the formation of biomolecular condensates. This study reports on the characterisation of a conserved domain located in the C-terminal region of the androgen receptor N-terminal domain, at atomic level, using nuclear magnetic resonance spectroscopy. This proline rich region exhibits extensive conformational heterogeneity driven by highly populated cis proline conformers that are stabilised through interactions with adjacent aromatic residues. We demonstrate that the cis-proline population is modulated by phosphorylation as well as cancer-associated mutations. This suggests that proline driven conformational heterogeneity at the C-terminal region of androgen receptor N-terminal domain is involved in the regulatory function of this transcription factor.

Protein structures encode non-local contact organization, but static coordinates do not directly quantify how a contact network responds when effective stabilizing interactions are strengthened or weakened. Here we introduce Contact-Network Responsiveness (CNR), a structure-derived framework that converts residue-level protein coordinates into density-controlled and topology-corrected response descriptors. The method is structure-source agnostic and can be applied to experimentally determined PDB structures, AlphaFold models, or other predicted structures; here, human AlphaFold models serve as the high-coverage structural substrate. Across 22,167 valid human protein structures, hydrophobic non-local contact density defined a nearly exact Bethe mean-field baseline for the conformational susceptibility threshold. A graph-aware residue-level extension then revealed systematic topology-dependent deviations from this density-only prediction. We define a topology correction ratio, [Formula], which separates topology-facilitated, density-dominated and topology-suppressed contact-network response regimes. CNR descriptors were associated with curated DisProt disorder annotations and broad-coverage UniProt/MobiDB-lite disorder fractions, supporting the interpretation that CNR captures a structural organization axis related to non-local contact availability and responsiveness.

Human glutathione synthetase (hGS) is a negatively cooperative ATP-grasp enzyme that catalyzes the final step in the biosynthesis of glutathione, a tripeptide antioxidant critical for life. hGS functions as an obligate homodimer with one active site per subunit; the two active sites are separated by [~]40 Angstroms. How ligand binding in one subunit reshapes the distant partner active site has remained a central unresolved question in understanding hGS regulation. This study provides the first atomistic model of ligand-dependent inter-subunit communication underlying negative cooperativity in hGS. Using atomistic simulations and dynamical network analysis, this study reveals how reactant- and product-bound states remodel the empty partner active site, redistribute inter-subunit interactions, and organize long-range communication between the two active sites. The product-bound/partner-empty state displayed a larger and less hydrated empty active site, demonstrating that ligand identity in one subunit alters both the geometry and solvent environment of the opposite site. Changes in ligand-dependent interactions are distributed across the dimer interface, with prominent contributions from the 42-46 interface region, the 11-30 region, and the 212-236 helical/interface region. Suboptimal path analysis shows product- and reactant-bound states share a communication scaffold, with 64.1% of transmission residues common to both pathways, 30.8% product-specific, and 5.1% reactant-specific. Together, the present results establish a detailed structural framework for hGS negative cooperativity in which ligand binding remodels a distributed allosteric network linking substrate-binding loops, the dimer interface, and the partner active site. More broadly, this work demonstrates how atomistic simulations can resolve long-range active-site coupling in multimeric enzymes and provides a foundation for experimental tests of allosteric transmission in hGS.

1.Virus-like particles (VLPs) formed from the murine polyomavirus major capsid protein VP1 are widely used as vaccine antigens and are being explored as nucleic acid and drug delivery vehicles. However, the factors controlling distinct VP1 capsid morphologies remain unclear. We investigated in vitro VP1 assembly with tRNA across NaCl concentrations of 0.15-1.0 M and tRNA mass ratios of 1:1-1:80 (w/w) using SEC-HPLC, transmission electron microscopy, and dynamic light scattering. Two competing assembly pathways were identified. At low ionic strength ([≤]0.15 M NaCl), nucleic acid-templated assembly produced compact, tRNA-filled T=1 VLPs ([~]28-30 nm). Assembly was maximal at tRNA ratios of 1:10-1:20, whereas excess or insufficient tRNA reduced yields. Increasing NaCl to 0.30 M lowered T=1 yields by 68-97%, and no T=1 particles were detected at [≥]0.5 M NaCl. Conversely, high ionic strength ([≥]0.5 M NaCl) promoted template-independent formation of hollow T=7 VLPs ([~]55-60 nm). T=7 assembly was inhibited by tRNA and was highest without nucleic acid. At 1.0 M NaCl, reducing the tRNA ratio from 1:40 to 1:80 increased T=7 yield more than 11-fold, while removing tRNA produced the greatest assembly efficiency. Kinetic analyses further showed that VP1 concentration and ionic strength regulate nucleation and assembly rate in the template-driven pathway. These findings show that electrostatic interactions govern pathway selection between tRNA-templated T=1 assembly and salt-driven T=7 self-assembly, providing practical guidance for controlling capsid morphology and cargo loading in VLP-based applications.

Nuclear organization is dynamic and originates from the fundamental subunit of chromatin, the nucleosome. Post-translational modification of nucleosomal histones, particularly within intrinsically disordered histone tail regions, provides a dynamic regulatory mechanism of accessibility for chromatin-templated processes. While the epigenomic impacts of lysine acetylation and serine phosphorylation in the histone H3 tail are well-known, how these charge-altering post-translational modifications (PTMs) alter nucleosomal tail conformational dynamics remains incompletely characterized. Given that the functional implications of these PTMs are, at least in part, a consequence of modified nucleosome conformation, systematically cataloging the impact of histone PTMs on nucleosome dynamics provides crucial insight into both baseline cellular activity and epigenetic dysregulation that occurs in disease. Previously, our lab demonstrated that arginine citrullination mimetics lead to regional increases in H3 tail dynamics within nucleosome core particles. Here, we performed nuclear magnetic resonance spin relaxation experiments to investigate the effects of lysine acetylation and serine phosphorylation on H3 tail picosecond-nanosecond (ps-ns) dynamics. Using lysine-to-glutamine and serine-to-glutamate mutations as acetyllysine and phosphoserine mimetics, respectively, we found that these PTMs increase ps-ns conformational dynamics regionally around the PTM site, with a position-dependent effect. Additionally, we show that the type of PTM influences the extent of these increases: in general, the effect of mimetics trends in the order of phosphorylation [&le;] acetylation < citrullination, suggesting a tunable method for altering histone tail dynamics. Taken together, these results illustrate the role of nucleosome conformational dynamics in conveying the effects of epigenomic PTMs, elucidating a mechanism of the histone language.

HELB is a human helicase involved in DNA repair and replication that interacts physically with the single-stranded DNA binding protein RPA. ATP-dependent translocation of HELB along ssDNA results in the active displacement of RPA molecules and the formation of ssDNA loops, suggesting that HELB contains at least two DNA binding sites. In this work, we investigated the role of HELB-specific structural elements in facilitating interactions between HELB and RPA-coated DNA. We show that a predicted OB-fold in the N-terminal region of the protein is important both for loop extrusion and RPA displacement. We confirm that a HELB-specific-motif within the RecA-like helicase/translocase domains is critical for binding RPA in solution but that, once HELB is bound to ssDNA, is dispensable for RPA displacement. We propose a model for RPA displacement in which both structural elements play important roles in the recruitment and activation of HELB at RPA-ssDNA filaments. SIGNIFICANCE STATEMENTSingle-stranded DNA generated during replication and repair is rapidly coated by replication protein A (RPA), creating a protected filament that must nevertheless remain accessible to DNA-processing enzymes. We show that the human DNA helicase HELB uses two specialised structural elements to overcome this problem. A HELB-specific RPA-binding motif promotes recruitment, whereas a predicted OB domain enables DNA looping and efficient RPA displacement. Removing these elements modestly enhances activity on naked ssDNA while impairing loading and activation on RPA- coated DNA. This reveals how helicase accessory domains restrain inappropriate motor activity while targeting the intended nucleoprotein substrate. Human HELB variants linked to reproductive ageing map to these regulatory regions, suggesting a connection between impaired RPA dynamics and reproductive health.

Neuroendocrine cells communicate with other cells by releasing neurotransmitters or hormones by exocytosis, which involves SNARE-mediated fusion between secretory vesicles and the plasma membranes of the secreting cells. In neurons two plasma membrane SNARE proteins, Syntaxin-1a and SNAP25, join with the vesicle membrane SNARE protein Synaptobrevin-2 to form a four-helix bundle, which drives membrane fusion. The assembly of these SNAREs, which is highly orchestrated in cells, has been intensely studied in solution using fragments of the SNARE proteins without their transmembrane domains or lipid anchors. However, in cell and model membranes, Syntaxin and SNAP25 are known to oligomerize and cluster, and little is known about how clustering affects their incorporation into SNARE complexes. In cells, the SM protein Munc18 has been implicated in aiding secretory vesicle docking and facilitating SNARE complex assembly through its interactions with Syntaxin. To understand how Munc18 orchestrates SNARE complex assembly on membranes, we employed protein reconstitution in model membranes as well as biochemical and biophysical assays to show that lipid-dependent oligomerization of Syntaxin affects Munc18-Syntaxin binding and SNAP25 insertion into the plasma membrane acceptor SNARE complex. We showcase the consequences of the different modes of Munc18-Syntaxin and SNAP25 interaction on Syntaxins oligomerization and orientation relative to the membrane surface, as well as on docking and fusion of purified insulin granules. We also determined low-resolution structures by cryoEM in nanodiscs and on the surface of proteoliposomes of membrane-bound assembly states of Munc18/Syntaxin and Munc18/Syntaxin/SNAP25 complexes.

Phosphoarginine (pArg) modifications direct proteins for proteolytic destruction by ClpC1P1P2, an essential mycobacterial protease that has emerged as a promising antibacterial drug target against Mycobacterium tuberculosis. The broader regulatory landscape surrounding pArg is poorly understood. Here, we establish a mechanistic connection between pArg binding and the activity of ClpC2, a non-proteolytic transcriptional repressor with homology to the ClpC1 N-terminal domain. Biophysical studies reveal that ClpC2 forms concentration-dependent higher-order oligomers that bind cooperatively to operator sequences in the clpC2 promoter. A high-resolution crystal structure of the Streptomyces thermoviolaceus ClpC2 C-terminal domain reveals a conserved dimerization interface mediated by a C-terminal helix, which is sterically disrupted by pArg binding. Consequently, we find that binding of pArg, as well as some ClpC1-targeting antibiotics, disrupts ClpC2 oligomerization, dissociates ClpC2 from its operator DNA, and relieves transcriptional repression in vitro. Moreover, comparative analysis of clpC2 promoters with single versus dual operator sites predicts differences in regulatory sensitivity across mycobacterial species. Together, these findings establish ClpC2 as a pArg-responsive sensor capable of mechanistically linking elevated pArg levels to downstream transcriptional regulation.

Fic enzymes mediate diverse post-translational modifications, including adenosine monophosphate (AMP) transfer and removal, referred to as AMPylation and deAMPylation, respectively. We identified the prokaryotic translation elongation factor Tu (EF-Tu) as an AMPylation target of the Fic enzyme SoFic. SoFic can constitutively reverse EF-Tu modification via deAMPylation whereas AMPylation depends on SoFic homodimerization. The complex crystal structure between SoFic and EF-Tu confirms a conserved target binding mode across evolutionary distant Fic enzymes. AMPylation disrupts EF-Tu's regulatory switch-I region, causing translational inhibition. SoFic furthermore binds to its promotor DNA, suggesting a dual function as transcriptional and translational regulator in bacterial cells. Together, our structural and biochemical data provide valuable insights into the functional and regulatory diversity of Fic enzymes.

An emerging result in the relationship between amino acid sequence and protein solubility is a preference, on average, for lysine over arginine in more soluble proteins. The termini of helices are known to be prone to partial unfolding, often employing N- and C-cap amino acids to maintain stability. Hypothesising that lysine/arginine differences in relation to solubility may be evident at helical termini, their propensities and predicted charge interactions in helices were examined. There is enrichment of lysine over arginine at helical termini in AlphaFold models of Escherichia coli proteins, more so (on average) for the most soluble proteins. Similar effects are seen for the sum of charged amino acids at helical termini. Regions other than helical termini also show correlation of lysine composition, and overall charged amino acid composition, with solubility. These results suggest that protein design protocols could improve solubility through targeting lysine enrichment in regions such as helical termini, in addition to the more conventional consideration of helix capping interactions.

Inverse-folding models can rapidly generate protein sequences compatible with a supplied backbone, but unconstrained redesign is poorly suited to enzyme and genome-editor-associated domains, where catalytic, substrate-proximal, and conserved structural regions must remain protected. In this paper, we present EditorForge, a modular constraint-and-audit suite for editor-domain protein redesign that wraps fixed-backbone inverse folding with explicit design masks, fixed-position enforcement, active-site-proximity auditing, active-site-shielded regeneration, and downstream structural quality control. Using full-length Moloney murine leukemia virus reverse transcriptase structure 4MH8 (MMLV RT 4MH8) as a demonstration target, EditorForge first restricted redesign to a bounded 25-position envelope while fixing 428 residues. An initial audit detected active-site-proximal failure modes despite fixed-position integrity. Later, the Active Site Shield module then removed five unsafe design positions, replaced them with lower-contact alternatives, and regenerated candidates under stricter constraints. Post Shield Audit evaluated 24 regenerated candidates, all of which satisfied the hard sequence/mask and active-site-shield constraints. For the eight candidates that were selected or returned for structure-prediction/refolding quality control, Enhanced RefoldQC found that all 8 evaluated predicted structures passed the computational structure-QC screen. That said, the selected 8 candidates passed the computational structure-QC screen, with global C RMSD values of 1.2061-1.5555 {degrees}A, active-site C RMSD values of 0.4098-1.8397 {degrees}A, mutation-neighborhood C RMSD values of 1.3155-1.6848 {degrees}A, and average pLDDT-like confidence values of 94.87-95.11. In short, EditorForge provides a reproducible triage layer that converts general inverse-folding output into constrained and editor-specific candidate sets for downstream structural and biological review on top of existing structural prediction tools.

More than two thirds of proteins in the human proteome are predicted to contain intrinsically disordered regions (IDRs), which lack stable folded structure. IDRs are critical for biological regulation and organization, as targets for post-translational modifications, and as mediators of biomolecular condensates. To address the pressing need for better structural models enabling functional insight, we developed AlphaFlex to model fully atomistic conformer ensembles for proteins predicted to have IDRs, modeled in the context of AlphaFold folded domains and an implicit bilayer for transmembrane proteins. The AlphaFlex resource provides conformational ensembles of human proteins from the AlphaFold database with identified IDRs in the Protein Ensemble Database that is mirrored in UniProt. This transformative resource of AlphaFlex ensembles provides physically and biologically relevant full-length models for IDR proteins, including scaffold proteins, those with IDR:folded-domain interactions, regulatory and condensate proteins requiring exposed binding elements, conditionally folding IDRs, and transmembrane proteins containing IDRs.

Cowpox virus CPXV012 inhibits MHC class I antigen presentation by interfering with TAP-dependent peptide transport, but its membrane-dependent structural organization and dynamic behavior remain incompletely defined. Here, we investigated the conformational properties of CPXV012 in membrane-mimicking environments and in a model of the CPXV012-TAP complex. CPXV012 was divided into three peptide constructs corresponding to the N-terminal cytosolic region, transmembrane segment, and C-terminal ER-luminal domain. The peptides were analyzed by circular dichroism spectroscopy, multidimensional NMR spectroscopy, and molecular dynamics simulations, and the resulting structural information was integrated into a full-length CPXV012 model. CD spectra showed that CPX-E1 and CPX-C2 are predominantly disordered in aqueous solution but acquire ordered, mainly -helical features in DPC micelles. NMR analysis in DPC-d38 micelles provided residue-level assignments and structural restraints supporting restrained structure calculations for both peptides. In three independent 1 {micro}s molecular dynamics simulations of the CPXV012-TAP complex, CPXV012 preserved a reproducible two-helical organization. The N-terminal/transmembrane region behaved as a relatively stable structural element, whereas the ER-luminal segment showed greater local flexibility. Interface analysis indicated that CPXV012 contacts both TAP1 and TAP2, with recurrent interactions concentrated in the luminal Y47-I69 region and involving polar and charge-complementary contacts. These results support a model in which membrane-associated structuring positions CPXV012 for TAP recognition, while the flexible ER-luminal region forms the main TAP-interacting surface. This structural framework complements existing functional models of CPXV012-mediated TAP inhibition.

Protein-protein interaction (PPI) information is distributed across resources that differ in organism coverage, identifier systems, evidence models, confidence scores and access mechanisms, so assembling and comparing evidence for a protein requires source-specific queries, identifier conversion and extensive post-processing. We present KlinkPPI, a web server that retrieves, compares and exports PPI evidence from STRING, BioGRID, IntAct, CORUM, HuRI and Predictomes from a single query. KlinkPPI accepts UniProtKB accessions, NCBI Gene and Ensembl identifiers and gene names, and performs taxonomy-aware mapping to a common identifier space. Users can query individual proteins across all resources available for an organism, or retrieve organism-wide interaction sets. Results are presented per source so that database-specific evidence, annotations and confidence values are retained, while an integrated view exposes coverage and agreement between resources. KlinkPPI deliberately does not merge heterogeneous confidence scores, nor collapse functional associations, complex co-membership, binary interactions and structural predictions into a single consensus network. Results can be exported in PSI-MI TAB 2.8-compatible or Apache Parquet format with user-selected evidence fields. KlinkPPI is freely available at https://rappsilberlab.org/KlinkPPI/ and the source code at https://github.com/Rappsilber-Laboratory/KlinkPPI. Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=110 SRC="FIGDIR/small/742057v1_ufig1.gif" ALT="Figure 1"> View larger version (16K): org.highwire.dtl.DTLVardef@12dcaceorg.highwire.dtl.DTLVardef@15d9d0aorg.highwire.dtl.DTLVardef@f0fc39org.highwire.dtl.DTLVardef@13f172e_HPS_FORMAT_FIGEXP M_FIG C_FIG

Membrane contact sites are organized by protein assemblies that physically couple organelles and coordinate lipid metabolism, yet the structural principles that enable lipid exchange across these junctions remain poorly defined. At the nuclear-vacuolar junction (NVJ) in budding yeast, the tethering protein Mdm1 and its binding partner Nvj3 form a complex that regulates lipid metabolic pathways, but the structural features underlying their interaction have not been resolved. Here, we use AlphaFold-based complex prediction and comparative structural analysis to define the organization of Nvj3-Mdm1 complex assembly. We identify a high-confidence heterodimer in which conserved PXA and PXC domains generate an extended tunnel spanning both proteins. Tunnel analysis predicts a core hydrophobic conduit traversing the Nvj3-Mdm1 interface, consistent with a lipid-compatible architecture. Evolutionary conservation is enriched at the Nvj3-Mdm1 interface. The predicted conduit shares geometric and physicochemical properties with bridge-like lipid transfer proteins, including Atg2, Fmp27, and Hob2, suggesting that heteromeric tether assemblies may contribute directly to inter-organelle lipid transfer. Cophylogenetic analysis reveals coordinated coevolution of Nvj3 and Mdm1 across Saccharomycetes. Together, these findings define Nvj3 as a structural partner of Mdm1 and support a conduit-based model of lipid transfer at the NVJ.

The ParABS system orchestrates chromosome segregation in many bacterial species. The centromere-like parS sites serve as nucleation points for the initial binding of the ParB protein. Subsequent diffusion on adjacent, non-specific DNA regions (spreading) in the presence of CTP and binding of more ParB molecules along with DNA looping via ParB-ParB interactions bring distal parts of the chromosome into proximity. ParB interaction with the ParA-ATPase motor protein, then, drives genomic segregation. It has been shown that in some bacterial species, the ParB-parS complex undergoes phase separation into a condensate. However, the physico-chemical properties of such condensates and their response to forces, such as those they may face in the cell, have not yet been characterized. Performing turbidity measurements in the presence of CTP and various concentrations of DNA and physiologically relevant mono and divalent salt It was shown that Mg2+ facilitates, while K+ concentrations higher than [~]20 mM disfavors, condensate formation. Microrheology measurements showed that condensates of ParB and DNA including parS sites (ParB-parS DNA) in the presence of CTP, are viscoelastic with a viscosity at Troom of [~]5 Pa{middle dot}s and able to quickly respond to deformations with a network relaxation time of 0.1 s. Additionally, fluorescence combined with force spectroscopy showed that mechanical disruption of ParB-DNA condensates in the presence of CTP requires [~] 5-7.5 pN of tension in the DNA, which is lower than the force required to stall a molecular motor such as RNA polymerase, but higher than the force required for the relocation of chromosomes and plasmids during segregation. These results support the idea that ParB-parS condensates dynamically rearrange at the molecular level while maintaining the cohesion necessary to sustain the drag force of segregation without interfering with genomic transactions. This physical mechanism could be the basis for the critical role of ParB-parS condensates in organizing and partitioning bacterial chromosomes. Statement of SignificanceCondensates formed by liquid-liquid phase separation enable cellular compartmentalization without the energy-costly production of a membrane enclosure and have been implicated in a wide array of cellular functions. Here, we analyze their chemo-physical properties in relation to their activity in bacterial chromosome segregation using the C. crescentus ParB protein and DNA containing parS specific binding sites. We find that the cohesion of condensates of ParB parS-containing DNA confers the mechanical stability necessary to pull genetic material towards the cell pole without interfering with genomic functions. We propose a mechanism to explain the function of ParB-parS DNA in chromosomal segregation.

A central challenge in molecular biology is understanding how amino acid substitutions modulate various features of protein function and stability. To illuminate the complexities of this relationship, high-throughput (HTP) assays are increasingly used to assess site-saturating mutagenesis libraries. A common downstream analysis is to average the set of twenty outcomes at each amino acid position for comparison with structural and evolutionary features. Average values clearly identify positions that tolerate most substitutions (neutral positions) and positions where most substitutions abolish activity (toggle positions). However, average values conceal the existence of rheostat positions, where different amino acid substitutions sample a wide range of outcomes. To quantitatively identify rheostat positions, we previously developed a histogram-based analysis that we here expand by: (i) incorporating new position classes observed in experimental studies of rheostat positions; (ii) formalizing a hierarchy of class assignments; (iii) refining error-based identification of neutral positions; and (iv) statistically assessing the robustness of class assignments to changes in experimental and computational parameters. RheoScale 2.0 is implemented in Excel and newly implemented in Python for facile integration with existing HTP pipelines; all parameters are customizable. Example analyses are shown for three HTP datasets of the SARS-CoV-2 papain-like protease. Results illustrate two aspects that influence interpretation of HTP data: First, position assignments (and substitution outcomes) depend highly on the measured feature. Second, many protein positions play multiple roles in the sequence-structure-function relationship. The recognition of varied position roles will advance understanding of pathogen evolution, protein engineering, and variant interpretation for personalized medicine. SummaryRheoScale 2.0 improves how high-throughput mutational data are interpreted by identifying protein positions where amino acid substitutions act like biological dimmer switches. By enabling more nuanced assignment of position behavior, beyond neutral or deleterious outcomes, this analysis framework advances studies of sequence-structure-function relationships and has broad relevance for understanding protein evolution, engineering proteins with desired properties, and interpreting variants linked to human disease. SOFTWARE AVAILABILITYhttps://github.com/liskinsk/RheoScale-calculator